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Abstract 
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The  purpose  of  nuclear  structure  research  is  to  understand  the  features  of  the  nuclear  force  that 
determine  the  various  ways  a  nucleus  can  behave  upon  excitation.  Theoretical  model 
calculations  are  compared  to  experimental  data  in  order  to  understand  which  models  better 
predict  different  nuclear  properties.  The  tellurium- 130  nucleus  has  two  valence  protons  with 
respect  to  a  closed  proton  shell  and  several  different  types  of  nuclear  structure  behavior  are 
thought  to  be  active. 

Experimental  data  were  collected  at  the  University  of  Kentucky  Nuclear  Structure  Laboratory 
using  a  technique  called  inelastic  neutron  scattering.  By  scattering  neutrons  off  of  130Te,  the 
nucleus  was  excited,  and  the  resulting  de-excitation  gamma  rays  recorded  as  Te  excited  states 
relaxed.  Through  analytical  techniques,  the  energy  level  scheme  has  been  constructed,  and 
spectroscopic  infonnation  such  as  lifetimes,  level  spins,  and  branching  ratios  has  been  obtained. 
The  behavior  of  the  130Te  nucleus  was  examined  from  the  viewpoints  of  the  General  Collective 
Model  and  the  Particle-Core  Vibration  Model.  The  General  Collective  Model  was  found  to  be 
inadequate.  The  Particle-Core  Vibration  Model  shows  promise,  but  the  best  technique  to 
improve  the  understanding  of  the  observed  level  scheme  and  properties  may  still  have  to  await 
the  availability  of  large-scale  Shell  Model  calculations. 
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I.  INTRODUCTION 


A.  A  Nuclear  Primer 

Many  people  have  taken  introductory  chemistry  courses  and  have  some  basic  knowledge  of  the 
atom  and  molecular  structure.  Few  have  taken  nuclear  physics  courses.  This  section  introduces 
the  reader  to  the  world  of  nuclear  physics  through  analogy  with  atomic  physics. 

The  atom  can  be  considered  a  collection  of  electrons  orbiting  in  a  force  field  created  by  a 
miniscule  point-like  central  object,  the  nucleus.  The  positively  charged  nucleus  attracts  the 
negatively  charged  electrons  through  the  “fundamental”  electromagnetic  force.  Because  the 
electrons  are  “bound”  to  the  atom,  the  system  is  quantized  and  the  electrons  are  found  only  with 
specific  energies  corresponding  to  specific  orbits.  If  energy  is  put  into  the  system,  electrons  can 
be  “excited”  to  higher  orbits  with  higher  energies.  In  the  hydrogen  atom,  where  there  is  only  one 
electron,  photons  are  emitted  (called  infrared,  optical,  and  ultraviolet  photons)  as  the  excited 
electron  falls  back  down  to  lower  states. 

In  a  multi-electron  atom,  many  electrons  can  be  involved  in  the  excitations.  The  whole  system 
again  can  only  exist  with  specific  energies  in  the  “ground  state”  or  in  other  specific  “excited 
states.”  Specifying  the  electronic  configuration  of  the  system  in  an  excited  state  becomes  much 
more  complicated,  and  uncovering  the  underlying  structure  of  an  excited  state  from  experimental 
observations  can  be  quite  a  chore  in  the  multi-electron  atom.  One  must  perform  a  model 
calculation  with  computer  codes  and  compare  its  predictions  with  measured  quantities. 

If  one  zooms  in  on  the  nucleus  of  an  atom,  one  will  find  that  it  is  not  a  point  object,  but  instead  a 
rather  large  blob  of  nucleons  (protons  and  neutrons),  with  a  radius  of  ~6  femtometers  for 
tellurium- 130  (130Te).  Each  nucleon  moves  in  the  attractive  force  field  created  by  all  the  other 
nucleons  through  residual  effects  of  the  “fundamental”  “strong”  force.  Much  like  in  the  atomic 
system,  because  the  nucleon  is  confined  to  a  region  of  space,  the  nuclear  system  is  quantized  and 
the  nucleons  occupy  specific  orbits.  If  energy  is  put  into  the  system,  the  nucleons  can  be  excited 
to  higher  orbits.  When  nucleons  fall  back  down  to  their  ground  states,  photons  (called  gamma 
(y)  rays)  are  emitted.  In  nuclei  with  high  atomic  numbers,  many  nucleons  are  involved  in  the 
excitations.  Still,  according  to  quantum  mechanics,  the  system  can  only  exist  in  certain  “excited 
states.”  Specifying  the  nucleon  configuration  of  the  system  in  a  specific  excited  state  becomes 
extremely  complicated,  and  uncovering  the  underlying  structure  of  an  excited  state  from 
experimental  observations  can  be  quite  difficult.  Again,  one  must  perform  a  model  calculation 
with  computer  codes  and  compare  its  predictions  with  measured  quantities. 

In  this  project,  we  excite  the  tellurium-130  nucleus  by  striking  it  with  a  fast  moving  neutron.  We 
measure  the  emitted  de-excitation  gamma  rays  and  detennine  each  excited  state’s  excitation 
energy,  total  angular  momentum  (a.k.a.  spin),  and  lifetime.  We  use  several  different  model 
calculations  in  an  attempt  to  uncover  the  underlying  nuclear  structure. 
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We  complete  the  introduction  with  a  slightly  more  detailed  introduction  to  the  project.  A 
description  of  the  accelerator  and  experimental  setup  is  provided  in  Section  II.  Section  III 
presents  the  specific  measurements  performed  and  the  analysis  techniques  that  allow  one  to 
extract  the  spectroscopic  information  from  the  data.  Section  IV  contains  experiment  results  and 
detailed  arguments  that  were  used  in  deciding  the  properties  of  each  excited  state.  The 
discussion  and  model  interpretations  are  presented  in  section  V.  Finally,  we  summarize  the 
results. 

B.  Nuclear  Structure 

The  role  of  modern  nuclear  physics  is  to  understand  how  the  interactions  of  dozens  and  dozens 
of  neutrons  and  protons  detennine  the  behavior  of  the  atomic  nucleus.  Just  as  the  van  der  Waals 
forces  between  molecules  are  the  result  of  residual  electromagnetic  interactions,  the  “nuclear 
force”  is  simply  the  result  of  residual  long-range  quark  interactions.  As  the  nucleus  is  a 
complicated  system  of  many  objects,  many  models  have  been  developed  to  try  to  understand  and 
predict  its  behavior,  and  this  research  continues  to  develop  the  existing  models  and  improve 
them. 

For  the  purposes  of  the  Introduction,  we  consider  three  classifications  of  nuclear  models:  the 
independent-particle  models,  the  collective  models,  and  a  combination  of  these,  called  hybrid 
models. 

The  independent-particle  models  [Kr88]  are  nearly  identical  to  the  quantum  mechanical 
treatment  of  atomic  (electron)  structure.  Nucleons  (neutrons  and  protons)  occupy  discrete 
orbitals,  and  may  jump  to  a  higher  orbital  if  they  absorb  energy  or  return  to  a  lower  orbital  by 
releasing  energy  via  gamma  rays.  Independent-particle  models  seek  to  explain  how  the 
individual  excited  states  are  constructed  from  nucleons  occupying  the  particular  orbitals  and  how 
the  combinations  of  orbitals  detennine  the  properties  (energy,  angular  momentum,  and  lifetime) 
of  the  excited  states.  However,  a  full  calculation  is  not  possible  in  the  case  of  the  nucleus 
because  a)  there  are  hundreds  of  particles  involved  and  b)  the  interactions  between  neighboring 
nucleons  are  very  strong.  All  independent-particle  models  have  a  theoretical  structure  somewhat 
like  that  of  the  larger  atoms  in  that  several  “valence”  nucleons  are  said  to  orbit  a  stable  inner 
core. 

The  collective  models  [Kr88]  describe  the  nucleus  as  a  fluid  of  neutrons  and  protons.  These 
models  seek  to  describe  oscillations  of  the  fluid  in  terms  of  the  normal  modes  of  oscillation  in  a 
liquid  drop.  Each  vibration  mode  has  a  definite  energy  associated  with  it,  and  this  gives  rise  to  a 
sequence  of  excited  states.  Rotational  excitations  also  occur  when  the  fluid  acquires  some 
deformation  and  thus  is  able  to  spin  like  a  top.  Again  the  models  seek  to  explain  the  properties 
(energy,  angular  momentum,  and  lifetime)  of  these  states  by  choosing  an  appropriate  potential 
energy  function. 

The  independent-particle  models  are  not  fit  to  describe  situations  where  large  numbers  (greater 
than  three  or  four)  of  nucleons  are  involved  and  the  collective  models  are  not  fit  to  describe 
situations  where  only  a  few  valence  nucleons  are  involved.  Hybrid  models  [Kr88]  seek  to 
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combine  features  of  both.  The  valence  nucleons  are  considered  to  orbit  a  core  that  can  undergo 
vibrations  and  rotations. 

Tellurium  was  chosen  for  this  and  previous  experiments  because  it  has  a  wide  range  of  abundant 
stable  isotopes.  Very  little  is  known  about  the  tellurium-130  isotope  nucleus  [Fi98,  BNL02, 
Ki02].  Although  it  is  naturally  abundant,  it  is  not  well  studied,  likely  because  it  is  not  accessible 
by  as  many  reactions  as  other  nuclei.  The  tellurium  nuclei  have  two  valence  protons  with  respect 
to  the  Z  =  50  (atomic  number,  number  of  protons)  closed  shell  and  a  range  of  neutron  numbers. 
Three  different  models  of  structure  are  thought  to  be  active  in  these  nuclei:  collective,  two 
valence  particle,  and  particle-hole  excitations  known  as  intruders.  Because  there  are  seven  stable 
Te  nuclei  each  with  an  even  number  of  neutrons  and  an  even  number  of  protons  (“even-even 
nuclei”),  one  can  study  the  evolution  of  these  excitation  modes  over  a  wide  range  in  neutron 
number,  allowing  for  good  quantitative  analysis  of  the  nuclei. 

The  ultimate  goal  in  nuclear  structure  research  is  to  understand  the  features  of  the  nuclear  force 
that  detennine  the  balance  between  the  various  collective-  and  particle-  like  behaviors.  This  is 
done  through  comparing  the  nuclear  models  to  data  taken  from  experiments  and  seeing  which 
models  work  better  under  which  circumstances.  For  this  Trident  project,  data  have  been  taken 
using  the  inelastic  neutron  scattering  process  on  130Te  and  recording  the  gamma  rays  that  are 
emitted  from  the  resulting  excited  nucleus.  Through  various  analytical  techniques,  a  scheme  of 
the  excited  state  energy  levels  has  been  constructed  complete  with  information  on  energies, 
lifetimes,  spins,  parities,  and  decay  branching  ratios.  This  scheme  was  then  compared  with 
model  calculations  to  see  how  the  different  models  describe  the  Te  nucleus.  Emphasis  centers 
on  understanding  the  interplay  between  particle  and  collective  features  and  on  the  aspects  of  the 
nuclear  forces  and  shell  model  orbitals  that  determine  the  relative  importance  of  each  model. 
Often  level  energies  are  more  important  for  evaluating  whether  a  model  has  sufficient  built-in 
complexity  than  distinguishing  between  physical  descriptions,  and  selection  rules,  lifetimes,  and 
transition  rate  infonnation  are  crucial  to  revealing  the  amplitudes  of  particle  and  collective 
components  in  the  wave  functions.  Results  from  these  Te  studies  will  be  combined  with 
previous  infonnation  on  the  other  lighter  tellurium  nuclei  120Te  [VaOl,  Va02],  122Te  [Co],  l24Te 
[Et97,  Wa98],  126Te  [Ta99,  TaOl],  and  128Te  [Ch97a,  Ch97b]  to  eventually  study  the  evolution  of 
these  structures  across  an  isotopic  chain. 
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II.  EQUIPMENT  AND  EXPERIMENTAL  PROCEDURES 

Inelastic  neutron  scattering  (INS;  n,n'y)  is  an  excellent  method  to  obtain  spectroscopic 
information  for  nuclear  excited  states.  Using  neutron  collisions  to  excite  nuclei  is  non-selective 
in  the  states  that  can  be  excited,  as  only  the  neutron  energy  and  the  angular  momentum  limit  the 
possible  states  that  can  be  activated.  The  high  quality  (n,n'y)  facility  at  the  University  of 
Kentucky  undergoes  continuous  improvement  (www.pa.uky.edu/~marcus/nukes.html).  Recent 
upgrades  have  made  operation  much  more  stable  and  user-friendly.  Special  detector  shielding 
and  gating  techniques  are  available  to  eliminate  the  large  neutron  backgrounds  present  in  neutron 
scattering  experiments. 

The  accelerator  at  the  University  of  Kentucky  is  a  7  MV  single-ended  Van  de  Graaff  with  a 
pulsing/bunching  system  located  inside  the  high-voltage  dome.  Proton  beams  emerging  from  the 
radio  frequency  source  are  pulsed  by  sweeping  across  a  chopping  aperture.  This  10-20  ns  wide 
chopped  pulse  immediately  passes  through  a  klystron  buncher  that  reduces  the  pulse  width  to  1 
ns  on  target.  The  repetition  rate  of  the  pulsed  beam  is  1.875  MHz.  Protons  are  accelerated  and 
impinge  on  a  gas  cell  at  the  end  of  the  beamline  containing  tritium  at  atmospheric  pressure  to 
produce  a  secondary  beam  of  neutrons  with  the  p  +  3H->  n  +  3He  reaction.  Typical  proton  beam 
currents  are  2nA  or  ~10  particles/second.  The  neutron  production  rate  is  approximately  1 0 
particles/second. 


Figure  1.  View  at  the  end  of  the  beamline.  Protons  arrive  through  the  pipe  to  the 
right  and  strike  the  tritium  cell  at  the  end  of  the  beampipe.  The  130Te  sample  is  in 
the  small  vertical  cylinder  in  the  center.  Layers  of  square  copper  plates  shield  the 
y-ray  detector  from  stray  neutrons  and  y  rays. 

The  sample  of  interest  (  Te)  is  hung  in  the  emerging  spray  of  neutrons  that  exits  the  gas  cell. 
For  our  experiments,  a  ~50  g  metallic  sample  of  >99%  enriched  lj0Te  was  borrowed  from  the 
Materials  Research  Collection  housed  at  Oak  Ridge  National  Lab.  The  University  of  Kentucky 
Nuclear  Lab  and  University  of  Dallas  were  responsible  for  arranging  the  sample  loan  and  sample 
rental.  Since  we  were  identifying  gamma-ray  transitions  for  the  first  time,  it  was  crucial  to  have 


a  sample  of  high  purity. 
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The  standard  detector  arrangement  employs  a  single  “>50%  efficiency”  high-purity  Germanium 
(HPGe)  detector  and  an  anti-Compton  shield.  The  Compton  shield  consists  of  a  Bismuth 
Germinate  (BGO)  crystal  with  six  photomultiplier  tubes  (PMTs)  which  surrounds  the  HPGe. 
Many  gamma  rays  do  not  deposit  their  full  energy  into  the  Germanium  crystal,  but  instead 
bounce  out  and  are  detected  by  the  surrounding  BGO  detector.  If  gamma  rays  are  detected  in  the 
BGO  and  HPGe  at  the  same  time  (in  the  time -resolution  of  the  detection  system),  then  they  are 
assumed  to  be  the  same  gamma  ray,  only  scattered.  Therefore  the  data  corresponding  to  that 
gamma  ray  event  is  meaningless  and  not  recorded.  The  detector  is  mounted  inside  graded 
copper,  lithium-loaded  paraffin,  and  lead  shielding  and  rests  with  its  cryostat  and  electronics  on  a 
detector  carriage  which  may  be  rotated  through  scattering  angles  up  to  150  degrees.  A  spectrum 
from  the  HPGe  detector  is  shown  in  Fig.  2. 
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Figure  2.  A  gamma-ray  spectrum  obtained  with  the  high-purity  Germanium  detector  - 
this  is  one  of  the  excitation  function  data  files  from  the  130Te(n,n’y)  data  runs  in  June. 

Approximately  150  lines  are  present  in  each  detector  spectrum. 

Data  is  recorded  by  use  of  a  Pentium  PC  based  -  Computer  Automated  Measurement  and 
Control  (CAMAC)  data  acquisition  system.  The  system  allows  both  online-sorting  and  event¬ 
mode  storage.  Pulses  from  the  HPGe  are  shaped  and  amplified  before  they  reach  the  analog-to- 
digital  converter  (ADC).  Based  on  the  recorded  energy  of  the  gamma  ray  captured  in  the  crystal, 
one  count  will  be  placed  somewhere  in  a  16384-channel  histogram.  The  channel  number-gamma 
ray  energy  relationship  was  calibrated  using  a  radioactive  radium-226  source  for  known  gamma 
ray  transition  energies.  The  detector  efficiency  (probability  that  a  y  ray  of  a  given  energy  stops 
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in  the  detector)  was  also  calculated  using  the  known  transitions  in  the  Ra  source.  Sufficient 
beam  time  was  assured  since  this  experiment  is  consistent  with  the  mainstream  experimental 
effort  of  the  nuclear  structure  group.  With  the  exception  of  several  typical  problems  with  the 
accelerator,  we  obtained  two  good  weeks  of  beam  time  during  the  summer. 

The  available  technical  support  in  the  physics  and  chemistry  departments  at  the  University  of 
Kentucky  is  excellent.  Services  include  a  full-time  accelerator  engineer  and  part-time  assistant. 
Fully  staffed  electronics,  machine,  and  glass  blowing  shops  are  available  free  of  labor  charges. 
The  accelerator  is  operated  by  the  users. 

Excitation  function  measurements  were  done  in  mid  June  2002.  While  in  Kentucky,  “accelerator 
watch”  was  stood  by  the  researchers.  This  duty  included  paying  attention  to  machine  operation, 
saving  data  periodically,  and  checking  levels  and  gauges  while  analyzing  gathered  data.  Angular 
distribution  measurements  were  done  during  the  first  week  of  August  2002.  At  this  time,  work 
began  on  the  excitation  function  analysis.  Advisors  Jeffrey  Vanhoy,  USNA,  and  Sally  Hicks, 
University  of  Dallas,  coordinated  the  use  of  the  accelerator  and  acquisition  system. 

Gamma-gamma  coincidence  measurements  done  on  natural  tellurium  by  Brian  Champine, 
USNA  ’98  in  1997  were  used  to  help  place  some  of  the  gamma  rays.  The  measurement  involves 
setting  up  multiple  detectors  around  the  target  in  close  proximity.  When  two  or  more  gamma 
rays  are  detected  in  close  temporal  proximity,  they  are  recorded.  This  can  show  which  gamma 
rays  are  part  of  the  same  decay  chain.  The  resulting  data  helped  to  confirm  several  suspicions 
about  the  level  schemes  while  making  placements  from  examination  of  the  excitation  functions. 


III.  ANALYSIS 
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A.  Peak-stripping  Codes 

After  the  data  have  been  taken,  the  process  of  analysis  begins.  The  choice  of  analysis  method  is 
important  because  some  methods  are  better  equipped  to  extract  accurate  and  precise  data  than 
others  under  various  circumstances.  The  analysis  package  needs  to  be  able  to  read  the  format 
from  the  University  of  Kentucky  data  collection  system,  plot  histograms  of  yield  vs.  ADC 
channel  number,  calibrate  the  channel  number/y-ray  energy  relation  using  known  radioactive 
sources,  fit  appropriately-shaped  peaks  to  the  histograms,  and  precisely  predict  the  shape 
characteristics  of  these  peaks. 

While  taking  data,  it  is  a  good  idea  to  analyze  what  you  have  in  order  to  determine  if  there  is 
something  wrong  with  the  accelerator  or  acquisition  system.  Several  programs  were  used 
initially  to  plot  data  “on  the  fly”  because  of  their  familiarity  and  therefore  ease  of  use.  These 
programs,  HYPERMET[Ph76]  and  PKS[Wa88],  are  Windows-based  programs  which  utilize  a 
user-created  input  file  and  the  data  file  to  fit  peaks  specified  by  the  user.  A  program  called 
SIMPLEFIT  was  also  used  to  integrate  the  area  of  the  histograms  given  user-defined 
backgrounds  and  peak  limits.  PKS  finds  its  own  backgrounds  and  looks  for  peaks  based  on  user- 
defined  window  limits  and  approximations  of  peak  centroids.  PKS  treats  Gaussian-shaped  peaks 
on  a  linear  or  quadratic  background.  Developed  originally  at  the  Naval  Research  Laboratory, 
HYPERMET  is  the  most  sophisticated  analysis  program  due  to  its  ability  to  model  peak  shapes 
but  does  not  have  an  extensive  graphical  interface  and  program  control  features.  For  the 
excitation  function  data  analysis,  a  program  package  commonly  used  at  the  University  of 
Kentucky  was  chosen.  The  package  is  a  set  of  Linux-based  programs  that  center  on  the  data 
stripping  program  “FITPIC”[UF96]  (see  Fig.  3). 


^  ■  •  mi  .if •  -  s  ■ _ ;  lv=  . -  ^ui 

Figure  3.  A  screen  shot  of  the  program  Fitpic  run  in  Linux. 
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FITPIC  has  relatively  good  capabilities  when  it  comes  to  fitting  the  peak  shapes  produced  by  the 
y-ray  detectors,  but  the  reason  it  was  chosen  above  the  others  is  because  its  graphical  control  is 
far  superior  to  that  of  PKS  or  HYPERMET.  This  makes  the  program  easier  to  learn  and  master, 
as  well  as  cuts  down  on  the  amount  of  time  spent  writing  input  files,  running  the  program, 
printing  results,  interpreting  difficulties,  re-writing  the  input  file,  and  repeating  the  process,  as  all 
of  this  can  be  done  at  the  same  time  when  using  FITPIC.  Another  advantage  of  using  FITPIC, 
with  its  graphical  control,  is  that  you  can  visually  examine,  interpolate,  and  change  multiple 
parameters  of  the  peak  shape  to  deal  with  variations  among  different  detectors  and  analog  to 
digital  converters.  The  satellite  programs  used  with  FITPIC  manage  detector  efficiency  and 
ADC  nonlinearity  corrections,  as  well  as  plot  the  data  in  a  form  necessary  for  proper  analysis. 

For  the  angular  distribution  and  lifetime  data  analysis,  the  Windows-based  program  PKS  was 
used.  After  using  FITPIC  initially  for  angular  data  stripping,  it  was  found  that  it  did  not  have  the 
ability  to  maintain  reliable  centroids  as  accurately  as  necessary  for  these  analyses.  After  peak 
areas  and  centroids  were  extracted,  spreadsheets  and  numerous  small  FORTRAN  codes  were 
used  to  further  reduce  the  data.  Angular  data-related  properties  (lifetimes,  branching  ratios, 
mixing  ratios,  etc.)  were  extracted  using  the  reduced  data  and  programs  dedicated  to  these 
individual  tasks. 


B.  Excitation  Functions 

The  most  logical  first  measurement  is  to  generate  the  set  of  excitation  functions.  An  excitation 
function  is  the  yield  of  each  gamma  ray  as  a  function  of  neutron  bombarding  energy.  These 
measurements  were  done  in  the  middle  of  June  2002.  To  get  the  set  of  neutron  energies,  the 
accelerator  was  run  at  each  energy  for  about  eight  hours  at  a  time,  and  the  beam  energy  was 
increased  by  approximately  90  keV  each  time.  This  was  done  from  beam  energies  of  1.89  to 
3.34  MeV,  recording  gamma  rays  with  the  detector  at  a  constant  90  degrees  to  the  beamline. 
Each  gamma  ray  that  is  produced  for  a  given  neutron  energy  has  a  yield  that  is  extracted  from  the 
area  of  its  peak  in  the  gamma  ray  spectrum.  When  the  yields  are  plotted  against  neutron  energy, 
the  minimum  neutron  energy  required  to  produce  that  gamma  ray,  the  “threshold,”  becomes 
evident.  Three  Te  gamma  ray  excitation  functions  are  shown  in  Fig.  4,  measured  in  beam 
energy  steps  of  ~90  keV,  enabling  a  y-ray’s  production  threshold  to  be  identified  to  <0. 1  MeV. 
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Figure  4.  A  selection  of  excitation  functions  from  the  30Te(n,n’y)  experiment. 
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Thresholds  are  determined  by  extrapolating  the  yield  plots  down  to  zero  yield,  and  seeing 
where  the  excitation  function  intersects  the  neutron  energy  axis.  This  is  the  energy  of  the 
bombarding  neutrons,  and  therefore  approximately  the  energy  of  the  y-ray’s  emitting  excited 
state.  From  the  excited  state  level  energies  obtained  through  analysis  of  the  excitation  functions, 
unidentified  gamma  rays  can  be  added/subtracted  to/from  these  level  energies  to  discover  new 
energy  levels.  For  example,  the  gamma  ray  with  energy  2403  keV  was  added  to  the  energy  of 
the  known  first  excited  state,  839  keV,  to  get  a  new  level  at  3242  keV,  which  coincides  with  the 
apparent  threshold  of  the  2403-keV  gamma  ray. 

Excitation  function  measurements  contribute  in  two  ways  to  the  level  scheme.  First,  knowing 
the  threshold  energies,  level  placements  are  made  for  all  the  apparent  transitions.  Second,  the 
shape  of  the  excitation  function  is  sensitive  to  both  the  spin  and  branching  ratios  of  the  parent 
state,  and  therefore  different  y-rays  with  the  same  excitation  function  shape  and  threshold  can  be 
identified  as  coming  out  of  the  same  state.  This  helps  by  aiding  in  placing  gamma  rays  whose 
excitation  functions  look  the  same.  The  complete  set  of  excitation  functions  is  given  in 
Appendix  A. 

C.  Angular  Distributions 

The  second  series  of  measurements,  taken  at  the  beginning  of  August  2002,  are  angular 
distributions.  Angular  distributions  are  measured  by  taking  spectral  data  at  a  series  of  detector 
angles  with  respect  to  the  beamline  at  a  fixed  incident  neutron  energy.  From  the  angular 
distribution,  two  pieces  of  experimental  infonnation  were  extracted  for  each  peak:  the  yield  of 
each  individual  y-ray  peak  versus  detection  angle,  and  the  centroid  energy  of  each  peak  versus 
detector  angle. 


The  angular  variation  of  an  individual  y-ray’s  yield  (Fig.  5)  will  tell  how  many  units  of  angular 
momenta  were  carried  away  from  the  nucleus  by  the  gamma-ray  photon.  This  infonnation  can 
be  used  to  extract  the  total  angular  momentum  (or  “spin”)  of  the  emitting  level  [Sh66].  The 
complete  set  of  angular  distributions  is  available  in  Appendices  B  and  C. 
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Figure  5.  Sample  angular  distributions  from  a  130Te(n,n’y)  experiment. 
The  line  is  a  fit  to  a  Legendre  polynomial  expansion. 


Level  spins  and  y-ray  mixing  ratios  are  determined  by  comparing  measured  y-ray  angular 
distributions  to  model  calculations.  Angular  distributions  are  fit  to  a  Legendre  polynomial 
expansion  generated  with  the  compound  nucleus  code  CINDY  [Sh66].  The  chi-squared  values 
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of  the  fits  of  Legendre  polynomials  to  the  angular  distribution  data  are  plotted  against  the 
arctangent  of  the  “mixing  ratio,”  or  how  much  of  each  type  and  multipolarity  (El,  Ml,  or  E2)  is 
present  in  a  measured  gamma-ray  transition.  The  multipolarity  is  how  much  angular  momentum 
is  carried  away  from  the  nucleus  by  the  detected  gamma  ray.  The  type  of  transition  is  either 
electric  or  magnetic  and  the  multipolarities  considered  are  1  and  2  [LiO  1  ] .  (The  arctangent 
function  is  used  since  the  mixing  ratio  can  go  from  negative  infinity  to  positive  infinity.)  Picking 
the  best  fit  (point  at  which  chi-squared  is  minimized)  from  these  plots  will  give  possible  spins  of 
the  source  excited  states  as  well  as  the  mixing  ratios.  Shown  in  Fig.  6  is  the  chi-squared  plot  of 
the  CINDY  comparisons  for  the  793-keV  gamma  ray. 


Engy  =  793 
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Figure  6.  The  value  of  chi-squared  for  fits  of  CINDY  model  angular  distributions  corresponding  to 
source  state  spins  of  1,  2,  3,  and  4  plotted  versus  the  tangent  of  the  mixing  ratio. 

Upon  inspection  of  the  figure,  one  can  see  that  the  value  of  chi-squared  is  minimized  for  the 
initial  level  which  produces  the  793-keV  gamma  ray  having  a  spin  of  4,  and  the  arctangent  of  the 
gamma-ray’s  mixing  ratio  is  near  -0. 15.  The  mixing  ratio  is  the  ratio  of  the  transition  amplitudes 
of  the  multipolarities.  Mathematically,  it  is  defined  as 
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5=<|L+1|>/<L>  [1] 

Where  <L>  is  the  reduced  transmission  amplitude  for  multipolarity  L. 

An  excited  state  may  decay  to  one  of  several  lower-lying  levels.  The  branching  ratio  (BR) 
indicates  what  fraction  of  time  a  particular  decay  branch  occurs.  Branching  ratios  were 
determined  by  comparing  the  angle-integrated  yields  of  all  the  gamma  rays  from  the  same 
excited  state.  Branching  ratios  of  low  energy  gamma  rays  must  be  corrected  for  a  process  known 
as  internal  conversion.  In  internal  conversion,  the  excitation  energy  is  transferred  to  the  atomic 
electrons  rather  than  being  released  in  the  form  of  a  gamma  ray.  This  process  is  generally 
negligible  for  transitions  greater  than  a  few  hundred  keV.  Internal  conversion  coefficients  were 
calculated  from  the  National  Nuclear  Data  Center  at  the  Brookhaven  National  Lab  website 
[BNL03], 


D.  Lifetimes 

The  angular  variation  of  an  individual  y-ray’s  energy  can  be  used  to  discover  the  lifetime  of  the 
emitting  level.  Lifetimes  were  extracted  with  the  Doppler-shift  attenuation  method  (DSAM) 
[Wi75,  Be94].  Although  there  are  other  techniques  for  detennining  lifetimes,  DSAM  offers  a 
good  approximation  for  lifetimes  between  one  femtosecond  and  one  picosecond.  Some  states  in 
130Te  are  known  to  have  lifetimes  longer  than  one  picosecond,  outside  the  bound  of  precise 
measurement  for  DSAM.  These  lifetimes  will  be  obtained  from  literature.  The  tellurium 
nucleus  recoils  when  it  is  struck  by  the  incident  neutron,  and  the  excited  residual  nucleus  emits  a 
y-ray  as  it  slows  down  in  the  target  material.  As  with  the  classical  Doppler  shift,  the  apparent 
frequency  of  the  photon  is  greater  in  front  of  the  moving  source  than  behind  it  (see  Fig.  7).  The 
complete  set  of  Doppler  shifts  is  available  in  Appendices  D  and  E. 

DOPPLER  SHIFTS 


Cos(fl)  Cos(0)  Co*(0) 

Figure  7.  Doppler-shift  data  for  three  transitions  in  130Te.  The  experimental  value  of  the 
Doppler-shift  attenuation  factor  F  is  determined  from  the  slope  of  the  best  fit  line. 


At  the  recoil  velocities  in  these  experiments,  the  energy  centroid  varies  with  the  angular 
dependence  indicated  in  the  following  equation: 

Ey(0)=EyO  [  1  +  Fp  cos(0)  ]  .  [2] 

The  product  Fp  can  be  thought  of  as  the  average  velocity  at  which  the  excited  nucleus  emits  the 
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gamma  ray.  The  symbol  [!  is  the  initial  velocity  of  the  recoiling  nucleus  (relative  to  the  speed 
of  light),  Ey°  is  the  true  y-ray  energy,  and  F  is  the  attenuation  factor  which  takes  into  account  the 
slowing  down  of  the  recoiling  nucleus.  The  nucleus  will  slow  down  in  the  material  due  to 
energy  loss  in  the  electron  cloud  and  collisions  with  other  nuclei  in  the  sample.  The  ratio  of  the 
speed  of  the  average  recoiling  nucleus  to  the  speed  of  light  for  the  3.3  MeV  data  set  was 
calculated  to  be  p  =  0.00064.  The  attenuation  factor  F  carries  the  dependence  on  the  level 
lifetime.  We  use  the  formalism  of  Winterbon  [Wi75]  to  relate  the  observed  attenuation  factors  to 
lifetimes  as  it  produces  more  reliable  results  for  a  greater  variety  of  target  compounds.  If  the 
slope  is  steep,  then  F  is  large,  the  nucleus  decays  relatively  soon  after  it  has  been  struck,  and  the 
lifetime  is  short.  If  the  slope  is  flat,  then  F  is  small,  the  nucleus  has  slowed  down  significantly 
before  decaying,  and  the  state  lifetime  is  long.  States  with  lifetimes  of  a  few  fs  to  1  ps  can  be 
determined  for  the  Te  nuclei  using  the  beam  energies  available  at  Kentucky.  Figure  8  shows 
the  relationship  between  the  attenuation  factor,  F,  and  the  lifetime,  t. 


Figure  8.  The  Winterbon  curve  relates  the  measured  attenuation  factor  and  the  source  state’s  lifetime. 

Much  effort  has  been  placed  on  improving  the  DSAM  method  at  Kentucky  during  the  past  10 
years.  This  effort  has  focused  on  reducing  and  monitoring  peak  shifts  (changes  in  the  building 
line  voltage,  HVAC  fluctuations,  and  amplifier  gain  drift)  in  the  y-ray  spectrum  and  improving 
the  understanding  of  stopping  powers  of  ions  recoiling  in  a  variety  of  target  compositions.  Gain 
shifts  have  been  greatly  reduced  by  special  power  isolation  and  special  cabling  for  detectors  and 
spectroscopy  amplifiers.  Close  track  is  kept  of  even  minute  gain  shifts  by  calibrating  detectors 
using  radioactive  Co,  Eu,  and  "  Ra  sources  before  and  after  measurements  at  each  detector 
angle. 


E.  Level  Scheme  Discussion 
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Once  the  values  for  lifetimes,  spins,  energies,  and  ratios  were  determined,  the  results  were 
compiled  into  a  single  spreadsheet  for  transition  rate  calculations  and  comparative  analysis.  In 
this  section  the  experiment  findings  are  expounded  upon.  Presented  first  are  several  discussions 
of  gamma  rays  that  caused  some  difficulty  in  placement.  Following  are  the  discussions  of  the 
level  scheme,  including  the  level  energies  and  characteristics.  The  compiled  data  table  is 
included  in  this  section  while  the  excitation  functions  are  found  in  appendix  A,  the  angular 
distributions  are  found  in  appendices  B  and  C,  and  the  Doppler  shifts  are  included  in  appendices 
D  and  E.  Figure  9  is  a  visual  representation  of  some  of  the  more  interesting  levels  in  the  level 
scheme  and  transitions  between  them. 


> 
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Figure  9.  An  abbreviated  version  of  the  130Te  level  scheme.  Values  at  the  right  are  level  energies  in  keV  and  state 
spins.  Blue  lines  represent  transitions  to  the  gound  state,  red  lines  represent  transitions  to  the  first  excited  state,  and 
green  lines  represent  other  transitions.  The  violet  lines  represent  the  transitions  of  interest  involving  the  2146  keV 

energy  level  described  below. 
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Difficulties  in  gamma  ray  placement  into  and  building  the  level  scheme  can  come  from  many 
different  problems.  Many  times,  the  spectrum  peaks  for  a  particular  excitation  function  have  low 
yields,  and  there  is  much  scatter  among  the  data  points.  Sometimes  the  excitation  energy  of  a 
state  is  too  high  to  get  an  adequate  number  of  data  points  for  an  extrapolation  to  find  the 
threshold.  Occasionally,  there  is  background  radiation  at  the  same  frequency  as  the  gamma  ray 
that  must  be  taken  into  account.  On  rare  occasions,  there  are  two  gamma  rays  from  the  nucleus 
with  nearly  the  same  energy.  All  of  these  things  cause  difficulty  when  trying  to  detennine  the 
threshold  of  the  gamma  rays. 

One  of  the  most  perplexing  problems  was  that  of  the  2146-keV  level.  In  the  excitation  functions, 
I  observed  no  gamma  rays  with  a  threshold  of  that  energy.  Therefore,  since  I  saw  no  gamma 
rays  coming  from  that  energy,  I  thought  there  must  not  be  an  energy  level  there.  The  problem 
was  that  several  excited  states  appeared  to  decay  into  a  state  at  2146  keV.  Specific  evidence 
included  a  gamma  ray  with  energy  258  keV  and  a  threshold  at  about  2400.  Upon  examining  the 
excitation  functions,  I  saw  that  there  was  a  gamma  ray  with  energy  331  keV  that  appeared  to 
have  the  same  threshold  and  general  shape  as  the  gamma  ray  with  energy  258  keV,  but  would 
not  fit  coming  from  the  same  level.  At  this  point,  I  consulted  previous  data  on  Te  and  found 
that  there  was  indeed  a  level  documented  at  2146  keV  [BNL02].  The  recorded  spin  of  the  state 
was  7”,  which  would  explain  why  I  could  not  observe  the  state  directly.  The  inelastic  neutron 
scattering  excitation  process  used  for  the  experiment  would  not  be  able  to  impart  seven  units  of 
angular  momentum  into  the  nucleus,  and  therefore  any  population  of  that  state  would  be  indirect, 
as  via  decay  of  higher  states  into  that  state.  The  gamma  ray  with  energy  331  keV  that  appeared 
to  come  from  a  level  at  2404  keV  actually  is  from  the  2146-keV  state.  The  reason  for  the 
delayed  threshold  is  that  the  2146-keV  state  could  not  be  populated  until  the  2404-keV  state 
which  fed  it  was  populated.  Therefore,  their  thresholds  appear  the  same.  And  since  the 
population  of  the  2146-keV  state  depended  on  the  decay  of  the  2404-keV  state  (via  the  258-keV 
gamma  ray),  the  gamma-ray  shapes  could  be  similar. 

The  285-keV  gamma  ray  is  responsible  for  difficulties  when  determining  another  level 
placement.  Its  position  is  a  transition  from  the  2171-keV  level  to  the  1886-keV  level.  It  appears 
that  the  gamma-ray  peak  could  be  a  “doublet,”  or  two  Te  gamma  rays  with  the  about  same 
energy.  This  is  observed  when  two  different  decays  with  two  sets  of  initial  and  final  states  have 
about  the  same  energy  difference.  The  excitation  has  a  severe  kink  near  2430  keV.  This  could 
indicate  a  doublet  with  gamma  ray  coming  from  a  level  at  that  energy,  but  unfortunately,  there  is 
background  radiation  at  that  gamma-ray  energy.  The  germanium  crystal  can  produce  a  285-keV 
gamma  ray  itself,  so  the  second  apparent  threshold  will  remain  speculation. 

Several  transitions  were  identified  with  help  from  the  gamma-gamma  coincidence  data.  For 
example,  the  334-keV  gamma  ray  has  a  threshold  around  2440  keV,  and  could  be  from  2435, 
2500,  or  2472.  But  the  gamma  ray  is  in  coincidence  with  a  468-keV  (2101  keV-1633  keV) 
gamma  ray,  so  I  picked  the  transition  between  levels  at  2435  keV  and  2101  keV  for  the  334-keV 
gamma  ray. 

Another  help  in  distinguishing  gamma  rays  for  level  placements  are  the  shapes  of  their  excitation 
functions.  Two  excitation  functions  are  said  to  have  the  same  shape  if  the  yields  of  all  the  points 
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in  one  excitation  function  nonnalized  to  a  single  point  in  it  are  the  same  as  the  yields  of  all  the 
respective  points  in  another  excitation  function  normalized  to  the  respective  single  point  in  that 
plot.  Excitation  functions’  distinguishing  shape  characteristics  include  the  slope  at  different 
points,  unique  kinks  in  the  line,  and  how  smooth  the  line  is.  The  shape  of  the  excitation  function 
of  a  gamma  ray  is  characteristic  of  the  branching  ratio  and  the  spin  energy  level  it  came  from. 
For  example,  the  385-keV  gamma  ray  could  have  been  from  any  of  three  levels  around  2700,  its 
excitation  function’s  shape  was  closest  to  other  gamma  rays  from  the  2789-keV  level.  The  491- 
keV  gamma  was  placed  from  2895  keV  because  it  looks  similar  to  the  288-keV  gamma  ray, 
which  comes  from  a  level  at  2895  keV.  Similarly,  the  521-keV  gamma  ray  looks  like  the  2117- 
keV  gamma  from  the  2956-keV  level,  so  it  was  placed  as  well  from  that  level.  The  575  keV  and 
614-keV  gamma  rays  apparently  both  come  from  2714  keV.  A  characteristic  kink  in  the  2138- 
keV  level  gamma  rays  indicated  that  the  103 1-keV  gamma  ray  was  one  of  them. 

Some  gamma  rays  have  baffling  characteristics  and  the  transitions  are  not  well  understood.  For 
example,  the  468.1  keV  (2101  keV  -  1633  keV)  is  a  very  strong  line  with  a  strong  excitation 
function  whose  yield  grows  continuously  with  a  kink  at  2850  keV.  This  type  of  “growing” 
behavior  is  neither  normal  nor  expected.  The  level  at  2101  keV  is  a  popular  choice  for  a  final 
state  for  gamma  decays,  which  might  have  something  to  do  with  the  busy  excitation  function. 
Unfortunately,  the  only  decay  from  this  state,  the  468-keV  gamma  ray,  is  too  low  in  energy  to 
determine  the  lifetime. 

Sometimes  gamma  rays  show  up  at  inconvenient  places.  For  example,  a  753.52-keV  gamma  ray 
has  a  strong  presence  in  the  excitation  function  data,  but  it  was  not  observed  in  the  angular 
distribution,  perhaps  due  to  the  duration  of  individual  runs.  It  turns  out  from  examining  the 
accumulated  data  on  the  other  tellurium  isotopes  that  it  is  from  either  12STe  or  126Te.  A  few 
gamma  rays  from  trace  elements  of  other  tellurium  isotopes  in  the  sample  are  always  present. 
Analyzing  the  data  and  previously  know  literature,  a  gamma  ray  at  636  keV  also  appears  to  be 
from  126Te,  and  the  743.3-keV  gamma  from  128Te  [BNF02].  Background  radiation  from  other 
elements  in  the  target  room  can  also  be  found  in  the  detector  spectra.  Copper,  iron,  and 
gennanium  gamma  rays  are  the  most  common. 

The  techniques  outlined  above  were  used  to  place  the  approximately  170  observed  transitions 
into  a  scheme  of  nearly  100  levels  lying  up  to  3.4  MeV  in  excitation  energy.  Spectroscopic 
information  about  each  state;  the  spin,  lifetime,  decay  branching  ratios,  and  multipole  mixing 
ratios,  have  been  extracted  for  many  of  these  states. 

Placing  levels  into  an  energy  scheme  using  excitation  functions  and  coincidence  data  was 
relatively  easy  compared  to  determining  lifetimes  and  other  angular-dependent  spectroscopic 
properties.  The  yield  of  a  particular  gamma  ray  is  much  less  sensitive  to  angular  dependence 
than  it  is  the  amount  of  energy  deposited  in  the  sample  by  the  neutrons.  Also,  with  the  very 
small  changes  in  energy  observed  for  Doppler  shifts,  it  is  often  difficult  to  extract  a  meaningful 
lifetime.  These  difficulties  make  the  calculation  of  actual  transition  rates  extremely  challenging. 

Angular  distributions  were  done  at  neutron  energies  of  both  2.2  MeV  and  3.3  MeV  in  order  to 
get  more  reliable  data  for  the  lower  energy  levels,  as  decays  from  higher  energies  tend  to  have  an 
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effect  on  the  observed  quantities  measured  for  the  lower  energy  states.  Running  with  lower 
energy  neutrons  does  have  significant  drawbacks,  though,  the  major  one  being  that  often  there 
are  not  enough  detector  counts  to  get  good  statistics.  Good  statistics  are  necessary  to  reduce  the 
scatter  among  the  data  points. 

Lifetimes  were  detennined  for  approximately  50  of  the  identified  excited  states  of  Te.  Many 
of  the  states  for  which  lifetimes  were  not  calculated  simply  did  not  have  adequate  statistics  in  the 
angular  data.  Often  what  was  observed  appeared  to  be  a  negative  slope,  which  is  a  physical 
impossibility.  The  DSAM  requires  identification  of  the  transition  energies  to  <0.1  keV  for  each 
angle.  Radioactive  Ra  sources  were  measured  between  each  tellurium  run  to  ensure  the  most 
accurate  energy  calibration  for  each  angle. 

From  the  plots  of  gamma-ray  yield  versus  detector  angle,  branching  and  multipole  mixing  ratios 
were  extracted  for  approximately  100  transitions.  Angular  momenta  were  assigned  to  most 
states,  with  new  spins  being  assigned  to  almost  50  energy  levels.  As  with  the  lifetime 
determinations,  analysis  was  complicated  by  scatter  in  the  data  points.  Spin  and  ratio  assignment 
relies  heavily  upon  the  quality  of  the  fit  of  CINDY’s  Legendre  polynomial  for  that  gamma  ray, 
as  well  as  the  fits  of  other  gamma  rays  from  the  same  state.  In  many  cases,  it  was  impossible  to 
determine  definitely  the  spin  of  a  state. 

Transitions  to  zero-spin  states  that  could  be  identified  as  having  a  bowed-up  or  bowed-down 
angular  distribution  identified  their  source  states  as  having  a  spin  of  two  or  one,  respectively,  as  a 
gamma  ray  can  only  carry  one  or  two  units  of  angular  momentum.  This  determination  did  not 
require  the  use  of  CINDY  as  these  transitions  are  purely  El  or  E2  in  type  and  multipolarity. 
Other  transitions  required  the  use  of  the  CINDY  code  because  the  different  possible 
combinations  of  angular  momentum  vectors  to  yield  the  difference  in  angular  momentum 
between  two  states  must  be  compared  to  see  which  is  most  plausible.  Decays  from  states  with 
zero  spin  have  flat  angular  distributions.  Therefore,  the  yields  of  the  1 125-keV  gamma  ray  from 
the  second  spin  zero  level  at  1964  keV  were  used  to  normalize  all  data  points. 

F.  Transition  Rates 

Reduced  transition  probabilities,  “transition  rates,”  for  different  multipole  transitions  were 
calculated  for  gamma  rays  from  the  spectroscopic  information  discussed  above.  Suitable 
references  for  this  material  are  Firestone  [Fi98]  and  Lilley  [LiO  1  ].  A  state’s  transition 
probability,  P(level),  is  defined  as  the  reciprocal  of  the  mean  lifetime,  1/x.  The  partial  gamma- 
ray  transition  probability  is  defined  by: 

Py(XL)=P(level)  N,(XL)  /  (X  Nd)  [3] 

Where  P(level)  is  the  transition  probability  for  the  entire  state,  Ny  is  the  intensity  of  the  gamma 
ray  of  interest,  X  Nd  is  the  sum  of  the  intensities  of  all  of  the  transitions  from  the  level.  The 
argument  XL  denotes  the  radiation  type,  X,  which  is  either  electric,  “E,”  or  magnetic,  “M,”  and 
the  multipolarity,  L,  which  is  1  or  2.  This  equation  can  more  simply  be  expressed  as  the 
branching  ratio  divided  by  the  lifetime,  BR/t.  For  pure  transitions  (El  or  E2),  this  equation  is 
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sufficient  for  the  partial  gamma-ray  transition  probability.  The  equations  for  a  single  decaying 
state’s  gamma-ray  transition  that  is  mixed-mode  (E2  and  Ml  present)  are  given  by: 

Py(M  1  )=P(level)/[  1  +S2+ay(M  1  )+82aY(E2)]  [4] 

Py(E2)=P(level)/[  1  +S'2+aY(E2)+S'2aY(M  1 )]  [5] 


Where  8  is  the  mixing  ratio,  and  aY  is  the  internal  conversion  coefficient,  which  for  essentially  all 
transitions  represents  a  minimal  effect.  If  more  than  one  gamma-ray  decays  from  a  state  that 
depopulates  via  a  mixed-mode  decay,  equations  [3]  and  [4]  or  [5]  must  be  combined.  The  result 
is: 


Py(Ml)=  BR/(t[1+S2]) 

[6] 

Py(E2)=  BR/  (t[  1+S'2]) 

[7] 

The  quantity  of  interest  is  the  reduced  transmission  probability  which  is  calculated  as: 


B{XL)  = 


L[(2L  +  l)!!]2/z 


r. 

nc 


8  tz(L  + 1) 


Ay 


Py(  XL) 


[8] 


where  the  units  for  energy,  length,  and  time  are  MeV,  fim,  and  second,  respectively.  The  B(XL) 
are  traditionally  converted  to  “Weisskopf  units”  in  order  to  easily  compare  the  transition 
probabilities  from  different  experimental  techniques.  The  conversion  factor  is  dependent  upon 
nuclear  mass  and  the  type  and  multipolarity  of  the  transition.  The  conversion  factors  are  as 
follows: 


B(E\)W=B(E1 )  /  1.65  [9. a] 

B(E2)W=B(E2)  /  39.1 1  [9.b] 

B(MI)W=B(M1)  /  1.79  [9.c] 

As  transition  rates  take  into  account  every  observed  property  of  the  gamma  rays,  they  are 
excellent  indicators  of  agreement  or  disagreement  with  model  calculations. 

All  spectroscopic  information  extracted  from  the  data  for  130Te  is  presented  in  Table  I. 

Transition  rates  are  presented  in  Weisskopf  units.  Uncertainties  are  presented  in  the  same  units 
as  their  respective  values  of  interest. 


Table  I.  Level  scheme  and  spectroscopic  information  for  the  excited  states  of  tellurium- 130.  Uncertainties  presented  have  been  symmetrized. 
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The  next  step  of  the  project  is  to  try  to  uncover  the  structural  features  of  the  Te  nucleus. 
Independent  particle  models  are  good  for  describing  systems  of  a  few  nucleons,  and  collective 
models  are  used  to  describe  large  nuclei  when  it  would  be  impossible  to  track  all  of  the 
individual  nucleons.  The  ultimate  goal  would  be  to  understand  the  features  of  the  nuclear  force 
that  detennine  the  balance  between  collective-  and  particle-  like  behavior.  To  look  at  the 
collective  features,  experimental  results  are  compared  to  the  calculations  of  a  collective  behavior 
model.  Likewise,  to  look  at  particle  features,  results  are  compared  to  particle  behavior  model 
calculations.  Hybrid  models  are  used  to  predict  relative  particle  aspects  of  larger  nuclei  by 
combining  aspects  of  both  collective  and  particle  models. 

Theorists  are  continually  developing  complex  codes  to  simulate  the  actions  of  excited  nuclei.  A 
number  of  these  treat  the  nucleus  in  the  “liquid  drop”  fashion,  proposing  that  all  excited  states 
can  be  reproduced  as  a  combination  of  vibrations,  deformations,  and  rotations  of  a  fluid  nucleus. 
An  older,  more  complete  traditional  fluid  model,  the  General  Collective  Model  (GCM)  has  re- 
emerged  [Gn69,  Gn71,  Tr91,  Va97,  Za97].  The  GCM  solves  the  problem  of  the  liquid  drop 
sloshing  around  in  a  deformed  potential  well.  In  the  GCM,  the  potential  energy  is  described  as  a 
function  of  two  defonnation  variables. 

The  independent-particle  models  are  not  fit  to  describe  situations  where  large  numbers  (greater 
than  three  or  four)  of  nucleons  are  involved  and  the  collective  models  are  not  capable  of 
describing  situations  where  only  a  few  nucleons  are  responsible  for  the  excitations.  Hybrid 
models  seek  to  combine  features  of  both.  The  Particle-Core  Coupling  Model  (PCM)  [He67] 
treats  the  tellurium  nuclei  as  two  valence  nucleons  orbiting  a  spherical  core.  The  two  valence 
nucleons  are  then  allowed  to  interact  with  each  other  and  with  the  core,  and  the  core  is  allowed 
to  undergo  simple  vibrations.  It  is  hoped  that  by  adjusting  the  strength  of  the  interaction  between 
the  valence  particles  and  the  core  (similar  to  tidal  forces),  one  can  “dial  up”  the  particular 
tellurium  isotope  of  interest.  Weaker  tidal  strengths  may  be  appropriate  for  the  vibrational-like 
120Te,  while  a  stronger  tidal  strengths  may  be  appropriate  for  130Te,  where  particle  characteristics 
are  more  pronounced. 

A.  Initial  Observations 

The  lighter  tellurium  nuclei,  ~  Te  and  '"Te,  have  groups  of  equally  spaced  levels  suggestive  of 
a  harmonic  oscillator.  The  heavier  nuclei,  126'130Te,  have  several  levels  suggestive  of  individual 
particle  orbits.  Looking  at  the  experimentally  developed  level  scheme,  it  appears  that  the  first  6+ 
state  is  too  depressed  to  be  part  of  a  harmonic  vibrator  scheme.  Also,  the  second  0+  state  is 
elevated  above  the  third  2+,  which  is  unusual.  Figure  10  depicts  the  evolution  of  a  few  excited 
states  across  the  isotopic  chain. 
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States  Across  the  Isotopic  Chain 
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Figure  10.  Selected  important  energy  levels  in  nine  tellurium  isotopes.  The  notation  of  the  legend  regards  the  first 

state  of  spin  2  as  “2.1”. 


Often  level  schemes  are  more  important  for  just  deciding  if  a  model  has  sufficient  built-in 
complexity  (since  the  Hamiltonian  parameters  are  adjusted  to  fit  the  observed  states)  rather  than 
providing  a  sensitive  test  of  the  wave  functions.  We  examine  the  electric  quadrupole  transition 
rate,  B(E2),  between  levels  to  evaluate  the  degree  of  agreement.  Formulae  for  determining  the 
transition  rate  from  the  experimental  lifetime,  multipole  mixing  ratio,  and  y-ray  transition  energy 
may  be  found  above,  as  well  as  in  other  references  [Kr88,  LiOl].  Comparisons  of  electric 
quadrupole  transition  rates,  B(E2),  for  several  transitions  of  interest  are  also  in  the  PCM 
discussion  section. 


B.  General  Collective  Model 
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The  General  Collective  Model  treats  the  nucleus  as  a  vibrating,  rotating  fluid  blob.  Rotations 
and  vibrations  are  quantum  mechanical  in  nature  and  correspond  to  quantized  quadrupole  and 
octupole  phonon  excitations.  The  radius  of  the  nucleus  is  expanded  by  the  following  function 
[He94]: 


R(0,cM)=Ro[1+£  a^(t)Y^(0,<|>)]  [10] 

where  a/p  are  expansion  coefficients  or  more  specifically  the  time-dependent  multipole 
deformation  parameters  corresponding  to  the  different  phonon  excitations,  and  are  the 
spherical  hannonics.  Typically  X  =  2  only.  Using  R(0,cf>,t),  kinetic  energy  terms,  elasticity 
tenns,  and  a  nuclear  potential  such  as  Woods-Saxon,  the  Hamiltonian  can  be  constructed  with 
independent  parameters  for  “inertia”  and  surface  stiffness. 

The  GCM  uses  what  is  known  as  the  Gneuss-Greiner  form  for  the  potential  energy  term  in  the 
Hamiltonian.  The  following  tenns  of  the  potential  were  used  [Tr91]: 

V(P,r)  =  C7^p2 cos(3rt  +  C4i/?4  [11] 

The  defonnation  variables  [!  and  y  determine  the  deformed  shape  of  the  nucleus  for  a  given  state. 
Changing  the  [!  and  y  coordinates  will  result  in  what  is  analogous  to  squashing  or  stretching  a 
balloon  in  the  z  direction  and  on  the  sides,  respectively.  The  GCM  code  takes  level  energies  and 
spins  given  by  the  user  and  searches  for  the  best  parameters  to  fit  those  inputs.  Then  the  user 
runs  the  code  in  a  different  mode  to  reproduce  a  level  scheme  with  these  “best”  parameters. 
The  quality  of  model  agreement  with  experimental  results  lets  the  user  determine  which  states 
exhibit  collectivity  and  which  states  do  not.  To  assure  that  the  minimization  routine  found  the 
absolute  minimum,  searches  were  made  starting  from  different  potential  parameters  selected 
from  other  tellurium  nuclei. 

When  the  GCM  code  analyzed  a  full  set  of  experimental  levels  to  determine  potential 
parameters,  the  values  it  determined  for  C2,  C3,  and  C4  were  1905.41,  20224.14,  and  4780.82 
MeV,  respectively.  There  are  other,  higher-order  tenns,  but  only  the  first  three  are  considered  to 
be  physically  justifiable.  The  GCM  was  also  asked  to  find  parameters  by  analyzing  a  reduced  set 
of  experimental  values;  the  “problem”  levels  were  taken  out.  The  potential  parameters  found 
with  less  search  levels  are  2042.498,  19721.71,  and  4786.108  MeV  for  C2,  C3,  and  C4, 
respectively.  Figure  1 1  is  a  comparison  of  the  experimental  level  scheme  and  the  two  GCM 
code  calculations. 
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Figure  11.  Comparison  of  the  experimental  130Te  level  scheme  with  GCM  calculations. 

Levels  are  arranged  in  columns  for  clarity. 

The  experimental  level  scheme  has  little  vibration-like  character.  The  position  of  the  first  state 
with  a  value  for  spin  of  2  and  a  positive  parity,  annotated  “2i+”,  was  correctly  placed,  but 
agreement  in  subsequent  levels  varied.  The  4j+,  22+,  and  02+  states  (1st,  2nd,  and  3rd  columns, 
respectively)  that  could  belong  to  a  two-quadrupole  phonon  multiplet  and  should  look  slightly 
degenerate  are  present  and  are  shifted  in  an  acceptable  manner,  but  the  fragmentation  of  the  6+, 
4+,  3+,  and  2+  states  of  the  three-phonon  multiplet  (second  column)  is  not  consistent  with  the 
patterns  of  anharmonic  vibrators  [Ca02].  The  three-phonon  states  are  intermingled  with  the 
two-phonon  states.  The  02+  is  shifted  above  the  2f  in  the  third  column.  The  6i+  state  in  the  first 
column  is  severely  depressed  from  that  of  an  equally-spaced  vibrator  scheme.  If  it  followed  a 
vibrator  scheme,  the  first  states  of  each  spin  would  be  evenly  spaced  (first  column),  as  additional 
identical  quadrupole  phonons  carry  the  same  energy.  The  sequence  of  states  that  would  be 
considered  a  part  of  the  emerging  y-  or  [1-  bands  (2nd  and  3rd  columns)  of  a  non-spherical 
vibrator  is  scrambled. 
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C.  Particle-Core  Coupling  Model 

We  compare  our  experimental  results  with  calculations  from  the  PCM,  which  treats  the  nucleus 
as  two  valence  protons  orbiting  an  inert  core  undergoing  harmonic  vibration  [He67].  The  PCM 
Hamiltonian  is  given  by  a  collective  term,  a  particle  term,  a  particle-particle  residual  interaction 
tenn,  and  a  proton-vibration  interaction  tenn.  The  ‘major  parameters’  for  the  PCM  model  are 
the  orbital  energies  for  the  valence  protons,  the  strength  of  the  pairing  interaction  between  the 
valence  protons,  the  energy  of  the  characteristic  core  vibration,  and  the  strength  of  the  interaction 
between  the  valence  protons  and  the  core.  Several  related  calculations  were  perfonned  some 
years  ago  on  various  tellurium  nuclei  by  Lopac  [Lo70]  and  Warr  [Wa98].  The  calculations  differ 
principally  in  their  choice  of  the  orbital  energies. 

Use  of  the  PCM  should  allow  us  a  better  understanding  of  the  particle  versus  collective  nature  of 
excited  state  structure  as  the  user  can  vary  the  interaction  and  strength  of  each  aspect.  Five 
single  particle  orbital  energies  can  be  defined  to  be  a  basis  for  the  particle-like  effects  in  the 
nucleus.  Multipole  phonon  energies  and  possible  combinations  may  be  varied  as  well.  The 
variable  parameters  2,2  and  C3  are  the  “coupling  constants”  of  the  quadrupole  and  octupole 
phonons,  respectively.  These  parameters  represent  the  particle-phonon  interaction  strength.  The 
variable  parameter  G  is  the  amount  of  interaction  between  the  particles,  the  “proton-proton 
pairing  strength.”  Although  the  older  PCM-generated  schemes  of  Lopac  and  Warr  were 
ultimately  chosen  to  represent  the  PCM  in  comparisons  with  experiment,  a  new  set  of  parameters 
were  also  independently  developed  in  an  attempt  to  reproduce  the  level  scheme  using  the  PCM. 

While  Te  has  two  valence  protons,  ~  Sb  only  has  one.  Since  it  has  an  extra  proton  beyond  the 
closed  shell  at  Z=50,  and  it  is  generally  more  difficult  to  break  a  closed  shell  than  it  is  to  excite  a 
single  valence  particle,  the  energy  levels  of  ~  Sb  probably  roughly  resemble  the  single-particle 
orbital  energies.  In  other  words,  the  low-lying  energy  levels  for  Sb  are  significantly  particle- 
like  in  nature.  Therefore,  the  tabulated  energies  of  the  excited  states  of  Sb  were  chosen  as  the 
five  single  particle  orbital  energies  as  follows: 

lg7/2:  0  MeV  2d5/2:  1.3  MeV  lhn/2:  1.13  MeV  2d3/2:  1.8  MeV  3si/2:  4.5  MeV 

The  “core”  nucleus  for  Te,  “  Sn,  has  a  closed  proton  shell,  and  its  first  excited  state,  the 
quadrupole  excitation,  should  exhibit  a  high  degree  of  collectivity.  On  such  a  basis,  128Sn’s  first 
excited  state  energy,  1.2  MeV,  was  chosen  as  the  130Te  quadrupole  phonon  energy.  Since  it  is 
thought  that  the  valence  protons  in  Te  are  involved  in  particle  behaviors,  only  the  neutrons 
were  allowed  to  take  part  in  the  quadruple  phonon  excitations,  and  no  octupole  phonon 
excitations  were  used  at  all.  Three  quadrupole  phonons  in  all  were  allowed. 

The  five  single-particle  energies  and  the  quadrupole  phonon  energy  were  entered  as  fixed 
parameters  in  the  PCM  input  files.  Often  the  parameters  C,  and  G  are  fixed  by  global  equations 
or  are  within  a  defined  range  dependent  upon  the  nuclear  number  and  mass.  For  this  experiment, 
the  energy  parameters  were  fixed,  and  these  coupling  strength  parameters  were  varied  in  order  to 
find  the  values  that  yielded  the  best  agreement  with  experiment  throughout  the  set  of  excited 
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states.  Seventy-two  sets  of  calculations  were  done  using  these  energies,  varying  the  pairing 
strength  parameter,  G,  from  0.04  to  1.00,  and  the  quadrupole  phonon  coupling  constant,  Ci,  from 
0.1  to  3.6.  Ratios  of  PCM-predicted  energies  to  the  experimentally  determined  energies  were 
plotted  as  the  dependent  variable  in  a  two-parameter  space.  Figure  12  illustrates. 


Figure  12.  The  ratio  of  the  PCM-predicted  value  over  the  experimental  value  is  plotted  as  a  function  of  both  the 
pairing  strength  and  the  quadrupole  phonon  coupling  constant.  The  boundary  between  the  red  and  blue  regions  is 

where  the  predicted  value  equals  the  experimental  value. 


Most  of  the  parameter  space  plots  look  similar,  and  it  is  easy  to  point  out  which  states  are  not 
easily  recreated  with  this  energy  set.  The  best  choice  appears  to  be  G=0.25  and  ^2=1.1.  These 
coordinates  were  closest  to  the  ratio=l  line  for  most  of  the  plots. 

From  the  plots,  one  observes  that  increasing  G  will  push  up  the  PCM-predicted  state  energies. 
Looking  at  the  plot  of  the  6i+  state,  one  can  see  that  the  optimal  value  for  G  is  less  than  what  was 
chosen,  and  based  on  the  results  that  the  state  was  calculated  to  be  too  high,  conclude  that  the 
proton-proton  interaction  is  less  than  predicted  for  this  particular  state.  Based  on  the  results  from 
the  GCM,  it  appears  that  the  6i+  state  is  indeed  particle-like.  Although  that  particular  instance  is 
an  extreme  case,  it  appears  that  many  of  the  levels  are  mostly  unaffected  by  a  change  in  the 
proton-phonon  interaction  strength.  Notable  exceptions  to  this  are  the  C>2+,  2\+,  and  4i+  levels, 
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which  exhibit  extreme  dependence  on  the  proton-phonon  coupling,  suggestive  of  mixed- 
character  states.  While  increasing  C2  for  those  states  reduces  their  energy,  the  opposite  appears  to 
occur  in  many  of  the  other  states,  including  the  later  2+  and  4+  states.  Figure  13  illustrates  the 
comparison  between  the  experimentally  detennined  level  energies  and  the  results  of  this  model. 


Figure  13.  Comparison  of  the  experimental  ,0Te  level  scheme  with  our  PCM  calculations. 

Levels  are  arranged  in  columns  for  clarity. 

Although  many  useful  conclusions  were  drawn  from  analyzing  the  level  energy  dependences  on 
G  and  £2,  the  resulting  level  scheme  for  the  chosen  values  does  not  reproduce  any  of  the 
irregularities  of  the  observed  tellurium  nucleus.  For  this  reason,  the  PCM  calculations 
previously  done  by  Lopac  and  Warr  were  chosen  for  further  evaluation  of  the  PCM.  Figure  14 
illustrates  the  comparison  between  the  experimentally  determined  level  energies  and  the  results 
of  the  model. 
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Figure  14.  Comparison  of  the  experimental  1  '°Te  level  scheme  with  the  Lopac  and  Warr  PCM  calculations. 

Levels  are  arranged  in  columns  for  clarity. 

The  proton  orbital  energies  in  the  Lopac  calculation  were  selected  to  best  reproduce  what  was 
known  about  the  tellurium  nuclei  in  the  year  -1970.  This  calculation  now  has  difficulty 
reproducing  the  current  level  sequences.  The  depressed  6+  state  is  not  reproduced.  The  Lopac 
calculation  is  more  representative  of  an  oscillator  with  a  y-deformation  dependent  potential 
(related  to  the  C3  term  in  the  GCM  potenial  of  equation  [11])  [Ca02]. 

The  proton  orbital  energies  in  the  Warr  calculation  were  based  upon  a  regional  shell  model  study 
by  Kisslinger  [Ki63]  and  the  Ph.D.  thesis  of  Wenes  [We83].  The  d5/2  and  d3/2  orbital  energies 
are  significantly  higher  than  those  in  the  Lopac  calculation.  The  resulting  level  sequence 
produced  does  a  reasonable  job  of  reproducing  the  ordering  of  states,  the  only  major  discrepancy 
is  the  elevated  position  of  the  second  2+  state  with  respect  to  the  4i+.  Because  the  PCM 
calculation  done  by  Warr  yielded  reasonably  sound  results  for  the  level  scheme,  it  is  known  the 
PCM  is  sufficiently  complex  to  model  the  Te,  and  further  evaluation  regarding  the  accuracy  of 
parameters  is  warranted.  The  electric  quadrupole  transition  rates  (B(E2))  are  therefore  compared 
(see  Table  II). 
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Table  II.  Selected  transition  rates  in  l30Te.  B(E2)  values  are  given  in  W.u.  Results  are  grouped  according  to  a 
pure  vibrational  model  interpretation.  Entries  denoted  are  not  measured.  Blank  entries  do  not 


apply- 


Initial 

State 

Final 

State 

Experiment 

Lopac [Lo70] 
Calculation 

Warr  [Wa98] 
Calculation 

Single-phonon  excitation 

2i+ 

o,+ 

15.1(3)  a 

31 

28 

Two-phonon  excitations 

4i+ 

2i+ 

- 

44 

34 

2T 

2i+ 

26.6(34) 

24 

35 

22+ 

o,+ 

0.016(2) 

1.7 

0.0015 

o2+ 

2i+ 

3.9(18) 

20 

20 

Three-phonon  excitations 

6i+ 

4,+ 

7.3(25)  a 

35 

21 

42+ 

4,+ 

13.0(38) 

6.9 

14 

42+ 

22+ 

not  observed 

13 

24 

42+ 

2i+ 

1.6(5) 

4.9 

7.5 

3i+ 

2i+ 

0.050(2) 

3.1 

0.36 

3i+ 

22+ 

15.4(76) 

32 

17 

3i+ 

4,+ 

22.7(112) 

4.7 

6.8 

23+ 

2i+ 

34.9(12)  or 
2.22(7) 

0.71 

3.5 

23+ 

22+ 

- 

0.78 

0.0001 

23+ 

4,+ 

- 

- 

- 

o3+ 

2i+ 

- 

0.0071 

0.12 

o3+ 

2i+ 

- 

9.7 

25 

Ground  state  decays  of  spin-2  states 

2i+ 

0i+ 

15.1(3) 

31 

28 

22+ 

o,+ 

0.016(2) 

1.7 

- 

23+ 

0i+ 

0.029(1) 

0.20 

- 

24+ 

0,+ 

0.37(3) 

0.0016 

- 

a)  Value  from  other  experiment  literature.  Ref  [VaOl]. 


There  does  not  appear  to  be  any  clear  correspondence  between  the  measured  rates  and  predicted 
rates,  either  in  absolute  size  or  in  a  simple  relative  large-versus-small  viewpoint.  The  actual  2i+  , 
02+,  and  6j+  states  are  not  as  collective  as  predicted  by  these  two  PCM  parameter  sets.  It  is 
difficult  to  make  structure  conclusions  for  the  remaining  levels.  As  an  inspection  of  the 
experimental  level  schemes  of  the  tellurium  nuclei  was  expected  to  be  dominated  by  a  mixture  of 
two-valence  proton  and  anharmonic  vibrational  structures,  more  exploration  with  parameter  sets 
in  the  PCM  is  warranted.  One  possible  approach  is  to  take  orbital  energies  from  experimental 
information  in  the  underlying  Sb  nuclei  again,  this  time  using  a  different  parameter  set  and 
phonon  energy[Co68,  Is67].  To  detennine  particle  effects  to  a  higher  degree  of  accuracy,  it 
would  be  wise  to  undertake  an  independent-particle  model  calculation.  Such  a  calculation  was 
not  started  for  this  research,  nor  has  it  been  attempted  by  more  experienced  shell-model  theorists. 


V.  SUMMARY 
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Nuclear  structure  research  is  not  a  trivial  process.  It  requires  much  time  and  effort  devoted  to 
planning,  experimentation,  and  analysis.  The  excited  states  of  tellurium- 130  have  been 
identified  from  gamma-ray  spectra  which  were  measured  following  inelastic  neutron  scattering. 
Data  were  collected  in  the  summer  of  2002  at  the  University  of  Kentucky  using  the  particle 
accelerator  and  data  acquisition  system  there.  Excitation  function  data  were  collected  in  June 
and  angular  distribution  data  were  collected  in  August.  The  data  have  been  analyzed,  and  the 
level  scheme  has  been  generated.  Over  half  (approximately  52%)  of  the  levels  were  previously 
undiscovered.  Level  and  transition  properties,  such  as  lifetimes,  spins,  mixing  ratios,  branching 
ratios  were  detennined.  Transition  rates  were  also  calculated.  Model  calculations  were 
performed  using  the  Particle-Core  Coupling  Model  and  the  General  Collective  Model. 
Comparisons  were  made  between  the  experimental  results  and  new  and  previously  existing 
model  calculations.  The  GCM  was  unable  to  reproduce  the  nuclear  structure  behavior  exhibited 
in  tellurium- 130.  A  previously  existing  PCM  calculation  was  found  to  reproduce  the  level 
scheme,  but  the  agreement  with  transition  rate  data  is  lacking.  As  the  experimental  level  scheme 
suggests  that  the  structure  is  dominated  by  a  mixture  of  two-valence  proton  and  anharmonic 
vibrational  structures,  a  more  thorough  exploration  with  new  parameter  sets  in  the  PCM  is 
warranted.  It  also  appears  that  more  particle-like  structure  is  active  in  several  levels  of  tellurium- 
130  than  was  previously  thought.  An  independent  particle  shell  model  calculation  would  be  the 
best  method  to  develop  a  better  understanding  of  the  nuclear  structure. 
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Appendix  A:  Excitation  Functions 
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The  excitation  functions  are  plotted  as  well  as  the  plot  of  gamma  ray  energy  versus  incident 
neutron  energy.  The  vertical  bars  on  each  point  represent  the  error  for  the  individual  point,  while 
the  horizonal  bars  in  the  energy  plot  represent  the  error  for  the  whole  set. 
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Gamma  energy  (keV)-  1766.04  Gamma  energy  {KeV)“  1756.04  Gamma  energy  (keV)-  176B.74  Gamma  energy  (KeV)*  1766.74 
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Gamma  energy  (koV)-  2117,42  Gamma  energy  (KoV)«  2117.42  Gamma  energy  (keV)—  2188.80  Gamma  energy  (KeV)-  21HB.B0 


76 


S'l  i  S‘0  0  8‘0  9'0  VO  CO  0  Z0‘0  SIO'O  l0-0j._0l.xs  0  CO  V 0  L‘0  0 

a«JD  psnjDuuaN  dud  psaj|DLiuofg  dud  pKi  olluon  mud  pea!|Duuo\j 


CO  SL'O  L*0  QO'O  0  Z  G'L  I  CO  □  L  G'Q  0  510  I'D  SO'D  0 

mjd  peejjBUJjafj  dojd  pes||Diiuoft  dsjd  p«j  dlujon  mjd  poh||dlu>3m 


BUZ  ^£1 LZ  £L  LZ  IBLZ  506 IZ  06  LZ  SCZZZ  CZZ Z  Q'TZZZ  VtZZ  S'Ct-CE  £*ZZ 

ffiueua  lxuujdc  duiludd  (A04)  X6jou3  owwog  (a^i)  cujujc^ 


Neutron  energy  (MeV)  Neutron  energy  (MeV)  Neutron  energy  (MeV}  Neutron  energy  (MeV) 


Gamma  energy  (koV)-  2270.5b  Gamma  energy  (KeV)«  2270.55  Gamma  energy  (keV)«  2280.26  Gamma  energy  (KoV)-  22B0.26 
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(A3)i)  ffiuoua  DVJUJD3  (/mbjQ  fbivw  dlulups  (A^>i)  XBueua  dwwds  (<\wi)  ABjous  dluiudo 
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Neutron  energy  (MeV)  Neutron  energy  (MeV)  Neutron  energy  (MeV)  Neutron  energy  (MeiV) 
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Appendix  B:  Angular  Distributions  at  3.4  MeV 


■UNCLE  (deg) 


ANCLE  (d#flj  ANGLE  (dm)  ANGLE 


Appendix  C:  Angular  Distributions  at  2.2  MeV 
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Appendix  D:  Doppler  Shifts  at  3.4  MeV 
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As  a  preface  to  the  collection  of  Doppler-shift  figures,  the  reader  is  reminded  that  the  energy 
shifts  involved  are  tiny.  In  the  centroid  equation 

Ey(0)=Ey°  [  1  +  Fp  cos(0)  ]  ,  [1] 

the  maximum  value  for  the  product  Fp  is  0.00064,  or,  we  are  looking  for  differences  in  energy 
«  0.06%.  The  HPGe  detector  resolution  is  2.2  keV.  As  a  result,  it  is  impossible  to  observe 
shifts  for  gamma  rays  with  energies  below  ~  900  keV  and  any  scatter  in  those  plots  is 
meaningless. 

Strong  peaks  in  the  spectrum  generally  provide  the  best  information  about  the  lifetime  of  the 
state.  For  weaker  peaks,  the  point  scatter  in  the  plot  is  a  result  of  poor  statistics.  Point  scatter  is 
also  large  when  there  are  two  close  peaks  in  the  spectrum. 
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Cos(th)  Cos(th)  Coa(th)  Coofth)  Cos(th} 
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Appendix  E:  Doppler  Shifts  at  2.2  MeV 
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